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ABSTRACT 

B1422+231 and PG1115+080 are gravitational lens systems produc¬ 
ing quadruple QSO images where there is real promise that time delays 
can constrain the Hubble constant. In addition, the lensing galaxies 
are both part of groups which can play an important role in modelling 
the lens potential. This article reports redshifts for the lensing galaxies 
and three neighboring galaxies in each of the two systems. B1422+231 
consists of a group ai z = 0.339 with a dispersion of 733kms“^, and 
PG1115+080 is a group at z = 0.311 with a dispersion of 326kms“^. 
One of the neighboring galaxies in the B1422+231 system turned out 
to be an emission line galaxy at z = 0.536, suggesting that QSO light 
passing through B1422+231 may have been subjected to lensing by a 
cluster at this more distant redshift. The velocity dispersion of the lens¬ 
ing galaxy in PG1115+080 is determined to be 281T25 km s“^ (1" square 
aperture), which is surprisingly large given the image splittings of 1.2" 
in that system. 

Subject headings: cosmology — distance scale — gravitational lensing — 
quasars: individual (Blj22+231, PG1115+080) 


^Based on observations at the W. M. Keck Observatory, which is operated jointly by the California Institute of 
Technology and the University of California 
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1. Introduction 

The differing paths followed by light around 
a gravitational lens leads to different time of 
flight. If the lensed source should vary, mea¬ 
surement of a time delay can transform di¬ 
mensionless redshifts into physical distances, 
hence provide a Hubble constant (Refsdal 
1964). To within factors of order unity the 
entire time of flight is simply the time delay 
divided by the square of the image splitting 
(in radians). In principle this is an extremely 
powerful method to learn about cosmology; in 
practice it has proven to be difficult to mea¬ 
sure time delays and model lenses accurately 
enough to challenge the accuracy of 10-20 
percent claimed by more traditional measures 
of Hq. Nevertheless, gravitational lens mea¬ 
surements of cosmology are critically impor¬ 
tant because they are potentially subject to 
far fewer biases than the usual distance lad¬ 
der. 

Time delays have now been measured for 
two lens systems: 0957-1-561 and PG1115-I-080. 
The difficulty in measuring a time delay and 
then rendering from it cosmological infor¬ 
mation is illustrated by 0957-1-561. A time 
delay was reported by Schild and Cholhn 
(1986), disputed by Press, Rybicki, and He¬ 
witt (1992), and hnally resolved in favor of 
Schild’s number by further observations by 
Kundic et ah (1996). Models of the system 
have also undergone improvement, with esti¬ 
mated values for Hq as low as 50 km s“^ Mpc“^ 
and as high as 90kms“^Mpc“^ (Grogin & 
Narayan 1996). 

PGl 115-1-080 has been observed by Schechter 
et al. (1997) who derived a time delay and es¬ 
timated a Hubble constant of 45 km s“^ Mpc“^ 
based on image and lens positions from Kris¬ 
tian et al. (1993) and a lens redshift of 


about 2 ; ~ 0.3 from Henry & Heasley (1986) 
and Angonin-Willaime et al. (1993). The 
observations reported here were recently an¬ 
ticipated by Kundic, Gohen, and Blandford 
(1997, KGB) who hnd 2 : = 0.311, and use this 
to derive Hq = 62 ±17. 

B1422-I-231 is a very well studied QSO and 
Ly-a system at a redshift of 3.62 discovered 
by Patnaik et al. (1992), with a lensing galaxy 
at a redshift of 2 ; = 0.647 according to Ham¬ 
mer et al. (1995). There is again a small 
group of three galaxies within a few arcsec- 
onds with unknown redshift. Although no 
time delay has been measured in this sys¬ 
tem, the photometry compiled by Keeton and 
Kochanek (1996) suggest that there may be 
enough variability in the QSO that it should 
be possible to do so. 

This article reports the hrst results in an 
ongoing program to measure the most ba¬ 
sic information necessary to exploit the grav¬ 
itational lenses: the redshift of the lensing 
galaxy, redshifts of galaxies clustered around 
the lens, and velocity dispersions of the stars 
within the galaxies. Although a small ingre¬ 
dient compared to measuring time delays and 
modelling potentials, it is an essential one. 

2. Observations and Reductions 

B1422-I-231 and PG1115-I-080 were observed 
on March 30 and March 31, 1997 using the 
Low Resolution Imaging Spectrograph (LRIS) 
(Oke et al. 1995) at the Keck H telescope on 
Manna Kea, along with calibration exposures, 
an observation of the cluster MS1358-I-62 to 
act as a velocity dispersion calibrator, and 
HD132737 and AGK2-I-14873 as radial ve¬ 
locity templates. The observations are sum¬ 
marized in Table 1. The sky was clear and 
the seeing was about 0.8" throughout both 
nights. Long slits of 1.0" and 0.7" were used 
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along with gratings of 300 and 600 l/mm, 
both blazed at SOOOA. The 300 l/mm grating 
was always rotated to provide coverage from 
3800-8700A, whereas the 600 l/mm grating 
was rotated either to cover 4800-7680A for 
the QSO observations or 3800-6080A for the 
templates. The spectral resolution was about 
7 . 9 A FWHM for the 300 l/mm grating, 4.65A 
FWHM for the 600 l/mm grating with a 1.0" 
slit, 3 . 59 A FWHM for the 600 l/mm grating 
with a 0.7" slit, and the scale along the slit 
was 0.211"/pixel. The template stars were 
guided smoothly across the slit in many lo¬ 
cations in order to have uniform illumination 
across the slit, to build up signal to noise, 
and to map out loci of constant slit position 
across the detector. The slit was rotated as 
illustrated in Figure 1, either to observe two 
companion galaxies simultaneously or else to 
cover the lensing galaxy while avoiding as 
much QSO light as possible. 

The spectra were reduced using software 
described in detail by Tonry (1984). The ba¬ 
sic steps are to flatten the images, remove 
cosmic rays, derive a wavelength solution as 
a function of both row and column using sky 
lines (wavelengths tabulated by Osterbrock et 
ah 1996), derive a slit position solution as a 
function of both row and column using the po¬ 
sitions of the template star images in the slit, 
rebin the entire image to coordinates of log 
wavelength and slit position, add images, and 
then sky subtract. A quadratic £t to patches 
of sky on either side of the object (including 
a patch between for the galaxy pair observa¬ 
tions) did a very good job of removing the sky 
lines from the spectra. 

As has been stressed by Kelson et al. (1997), 
measuring the velocity dispersion of a galaxy 
at a redshift significantly larger than zero 
must be done carefully, since the instrumental 


Fig. 1.— Illustration of the slit positions 
and galaxy idenihcations for B1422-I-231, 
PG1115-I-080, and MS1358-I-62. North is up 
and East is left for each diagram. 
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resolution of a spectrograph tends to a con¬ 
stant number of angstroms, whereas the red- 
shifted spectrum has been stretched. Hence 
a straightforward cross-correlation or Fourier 
quotient will underestimate the dispersion of 
the galaxy. Kelson et ah measured disper¬ 
sions in the cluster MS1358-I-62 using high 
resolution template spectra and very careful 
modelling of the spectrograph resolution de¬ 
rived by measuring sky line widths. The ob¬ 
servations here use the trick that the ratio 
between an 0.7" and 1.0" slit is slightly larger 
than {l + z) for these lenses, hence a template 
measured with the 0.7" slit will have almost 
exactly the same instrumental resolution as a 
galaxy at a redshift of 2 : ~ 0.3. This is borne 
out by the ratio of the measured spectral res¬ 
olutions: 4 . 65 A -=- 3 . 59 A = 1.30. This is ap¬ 
plicable only to the 600 1/mm observations, 
of course. 

As a test that this observing procedure will 
give correct dispersions, we also observed the 
cluster MS1358-I-62 for which Kelson et al. 
measure a dispersion of 305 ±5 for the central 
galaxy, Gl, and 213 ± 4 for the neighboring 
bright galaxy, G2 (Franx, private communi¬ 
cation). 

The analysis of the companion galaxies was 
straightforward since the spectra were of very 
high signal to noise and were uncontaminated 
by QSO light. In each case the spectrum 
was extracted, and cross-correlated with the 
template spectrum according to Tonry and 
Davis (1979) as well as being analyzed by 
the Fourier quotient method of Sargent et 
al. (1977). (The cross-correlation is more ro¬ 
bust in the case of low signal to noise, but 
at the signal levels here the two results are 
statistically the same, and are simply aver¬ 
aged.) For each spectrum the redshift, er¬ 
ror, and velocity dispersion were calculated. 


The redshift calculation included the entire 
spectrum, whereas the spectrum blueward of 
4000A (rest frame) was excised for the dis¬ 
persion calculation, since the calcium H and 
K lines are so broad that they are always 
problematic for dispersions. Dispersions are 
not shown for 300 1/mm observations, and 
they are not corrected for any aperture effects, 
hence correspond to an aperture of approxi¬ 
mately 1" square. 

Table 2 lists the redshifts, errors, velocity 
dispersions, errors, and cross-correlation sig- 
nihcance “r” values for each spectrum. 

2.1. B1422 

We were surprised to hnd an unmistakable 
emission line system adjacent to B1422-G2, 
offset by about 1.0" along the line towards 
B1422-G3. It showed very clear, extended 
emission lines of 3727 and 5007 in both ex¬ 
posures, and H/3 could also be faintly dis¬ 
cerned. Lacking any imaging information on 
this galaxy, we tabulate it as “B1422-Gx”. 

The lensing galaxy extractions require con¬ 
siderable care. In the case of B1422 we ex¬ 
ploited the fact that the A-B-G components 
of the QSO would put a more extended dis¬ 
tribution of light along the slit than the lens¬ 
ing galaxy. We therefore extracted three 
swaths of approximately 1"; the central one 
we dubbed “QSO-fgalaxy” and the two flank¬ 
ing swaths we added and dubbed “QSO”. Of 
course there is considerable galaxy light in 
the flanking swaths but the expectation was 
that the proportion of galaxy light would be 
less. We found a linear function of wavelength 
which, when multiplied by the “QSO” spec¬ 
trum, would match the Ly-a and G IV QSO 
lines in the “QSO-fgalaxy” spectrum. Sub¬ 
traction left behind a residue which was a 
clean enough lensing galaxy spectrum that we 
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could derive a solid cross-correlation redshift. 

In principle it would have been better to 
get a pure QSO spectrum for subtraction, 
and indeed we obtained a pure QSO spectrum 
by shifting the slit down on top of the A-C 
line. However, the QSO light in the compos¬ 
ite spectrum has a large gradient across the 
slit, hence a wavelength shift and large inten¬ 
sity variations. Time did not permit rotating 
the spectrograph 180 degrees and offsetting 
to the other side of the QSO to try to dupli¬ 
cate the QSO contribution to the composite 
spectrum. 

We see no emission in the B1422 lens- 
ing galaxy whatsoever, despite the claims by 
Hammer et ah (1995) that 3727 and 5007 
appeared at 6138A and 8247A. We excised 
the A and B bands at 7625A and 6885A 
from the spectrum and 130A centered on the 
Ly-a peak, and performed a cross-correlation 
against the template star. The cross-correlation 
came out with a correlation peak at 2 ; = 
0.3366 with an r valne of 4.2, which is signih- 
cant enongh that it is quite unlikely (probabil¬ 
ity around 1 percent) to be spurious. It is pos¬ 
sible to see H-l-K, Mg, and Na at that redshift. 
With the corroboration from the neighboring 
galaxies, we are qnite conhdent that this is 
indeed the redshift of the lensing galaxy. 

The spectra of the QSO, lensing galaxy, 
G3, G2, and Gx are plotted in Figure 2. 
The most prominent Frannhofer lines are la¬ 
beled for the lensing galaxy, and while they 
can plausibly be seen by eye, our trust in 
this redshift ultimately depends on the cross¬ 
correlation signihcance. 

2.2. PG1115 

As with B1422, we serendipitously picked 
np an emission line galaxy in the PG1115 sys¬ 
tem along the slit between G1 and G2, lying 


Rest wavelength (A) 



Fig. 2.— The spectra are shown for G3 (off¬ 
set by 500 for clarity), G2, Gx, QSO-l-lensing 
galaxy, “pnre” QSO, and extracted lensing 
galaxy (scaled by a factor of 10 for visibil¬ 
ity) in the B1422 system. These spectra have 
been Gaussian smoothed with a fwhm of 6A. 
The top axis shows rest wavelength at a red¬ 
shift of 0.34. The [OH] and [OHI] lines in Gx 
are labeled as well as the Fraunhofer lines in 
GL. 
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3.8" from G2. It displays emission lines of 
5007, 4959, H^, and H^, and we refer to it in 
Table 2 as “PG1115-Gx”. 

The slit for PG1115 ran over the top of 
the B QSO component, the lensing galaxy, 
and skirted the edge of the G component, cre¬ 
ating a bimodal distribution of light in the 
slit with the galaxy in the center. In this 
case the QSO light illuminated the slit uni¬ 
formly enough that the subsequent exposure 
centered on the A1-A2 components matched 
the QSO light from the B and G to a high 
degree of accuracy. The B and G peaks were 
separated by 1.8" at this position angle, and 
we extracted the inner 1" from between the 
two peaks and subtracted a scaled version of 
the A1-A2 spectrum to give us the lensing 
galaxy spectrum illustrated in Figure 3. 

Gomparison with the redshifts presented 
by Kundic et ah shows extremely good agree¬ 
ment: GL - 0.3100 KGB, 0.3098 here; G1 - 
0.3099 KGB; 0.3098 here; G2 - 0.3120 KGB; 
0.3123 here. As noted by KGB, these red- 
shifts are in good agreement with Henry & 
Heasley (1986), but are only marginally con¬ 
sistent with Angonin-Willaime at ah (1993). 

The spectra of the QSO, lensing galaxy, 
Gl, G2, and Gx are plotted in Figure 3. 

2.3. MS1358 

A pair of exposures of MS1358 identical 
to those of PG1115 yielded spectra of seven 
galaxies. The slit was aligned between the 
central galaxy, called Gl here, and the bright 
galaxy 19" to the SW (PA 204), called G2. 
The galaxies G3 and G4 were well centered 
in the slit, hence their velocities and disper¬ 
sions should be reliable. The three remaining 
galaxies were marginally caught by the slit 
and low signal-to-noise, so we do not report 
their velocities. Gomparison with Fabricant 


Rest ■wavelength (A) 

4000 4500 5000 5500 



Fig. 3.— The spectra are shown for Gl, G2, 
Gx, QSO-|-lensing galaxy, pure QSO, and ex¬ 
tracted lensing galaxy in the PG1115 system. 
These spectra have been Gaussian smoothed 
with a FWHM of gA. The top axis shows rest 
wavelength at a redshift of 0.31. Emission 
lines of 5007, 4959, and 4861 are apparent in 
the spectrum of Gx. 
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et al. (1991) for the three galaxies in com¬ 
mon (G1 = #140, G2 = #121, and G3 = 
#126) gives a mean offset of 5z = 0.0002 and 
a scatter of az = 0.0001, well within the error 
estimates. 

The comparison of velocity dispersions with 
Kelson et al. reveals no significant offset in 
the dispersions measured here — Gl: 305 ± 5 
versus 299 ± 22 here; G2: 213 ± 4 versus 
235±23 here. The mean of these ratios differs 
from unity by only 4%, indicating that the 
dispersions and errors reported here should 
be accurate. 

3. Discussion 

Gomputing means and standard deviations 
for the redshifts presented here, we hnd the 
B1422 lens is part of a group at Zg = 0.339 ± 
0.002 with a rest frame dispersion (i.e. stan¬ 
dard deviation of cz divided by 1 + z) of 
733kms“^ (3 galaxies). The PG1115 lens is 
part of a group at Zg = 0.311 ± 0.001 with a 
rest frame dispersion of 326kms“^ (5 galax¬ 
ies). Kundic et al. gave 165kms“^ but this is 
obviously a typographical error which should 
be 270kms“^, (conhrmed by Kundic, pri¬ 
vate communication). Using the four galax¬ 
ies observed in the MS1358 cluster we hnd 
Zc = 0.325 ± 0.001 with a rest frame disper¬ 
sion of 432kms“^ (4 galaxies). 

A few comments seem appropriate. The 
B1422 group appears to be quite massive, par¬ 
ticularly given the number of galaxies present. 
Deeper imaging would probably be worth¬ 
while. Despite the comparison Kundic et al. 
make between the PG1115 and Hickson poor 
groups (Hickson et al. 1992), it now looks as 
though the group is much more massive, and 
the dispersion of 326kms“^ is in quite good 
agreement with the models of Schechter et 
ah, who need a group dispersion of 383 kms“^ 


to provide the necessary shear to explain the 
image positions. Finally, the dispersion for 
MS1358 serves as a warning that small num¬ 
bers of galaxies near the centers of groups 
may not be a reliable measure of the over¬ 
all group mass. Fabricant et al. determine 
a dispersion of 1125 km s“^ for MS1358 from 
approximately 65 galaxies within 2 arcmin- 
utes of the BGG. Recomputing a dispersion 
from their redshifts but restricting to galaxies 
closer to the BGG than G2 we hnd a disper¬ 
sion of 625kms“^, demonstrating how a cool 
subsystem can exist in a massive cluster. 

The redshift found here for the B1422 sys¬ 
tem is signihcantly lower than that reported 
by Hammer et ah, so the lens and companions 
are not as highly luminous as had been sug¬ 
gested. This lower redshift is in quite good 
agreement with the redshift estimated pho¬ 
tometrically by Impey et al. (1996). The 
rather high cluster mass called for by Hogg 
and Blandford (1994) appears to be conhrmed 
here. The presence of the background galaxy 
at z = 0.536 at a distance of 9" from the lens 
may point to a background cluster which is 
providing signihcant multiple lensing. Very 
deep imaging of this source would help elu¬ 
cidate whether there is indeed a background 
cluster present. 

The dispersion for the lensing galaxy in 
PG1115 of 281kms“^ can be corrected ac¬ 
cording to the formula given by Jorgensen 
et al. 1995 to the standard metric aperture 
commonly used in fundamental plane studies: 
3.4" at the distance of Goma. This corrected 
dispersion is 293kms“^ and it seems surpris¬ 
ingly large. Perhaps the best way to see this 
is to consider the Einstein ring created by an 
isothermal lens of dispersion a aligned with 
the source. This ring will have radius rg given 
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by 


DdsDd 4:71(7"^ 

Ds ^ ^ 

Using redshifts of Zg = 1.722, Zd = 0.311, and 
f2o = 1) we find the first term is 394/i“^ Mpc, 
and the angular size distance of the lens 
is Dd = 580/i“^ Mpc, so that one arcsec- 
ond is 2.81/i“^ kpc. This then gives tq = 
A.72h~^ kpc = 1.68". However, the geomet¬ 
ric mean of the four QSO positions around 
the lensing galaxy is 1.17", which is the rea¬ 
son that the models of Schechter et ah prefer 
a dispersion for the lens of about 235kms“^. 
The position of the ring of images is quite 
robust to the details of shear and azimuthal 
image position, so it will take detailed mod¬ 
elling to understand how this high dispersion 
can £t in with the rest of the observables. It 
is conceivable that this dispersion is simply 
wrong, but the spectrum is quite high qual¬ 
ity, even given the necessity for QSO subtrac¬ 
tion, having a signal to noise of about 20 per 
angstrom. It is also possible that the disper¬ 
sion drops rapidly from the measured value of 
281kms“^ to something nearer 235kms“^ at 
the radius of 3.3h~^ kpc where the images lie. 
While not unheard of, this would be a sur¬ 
prisingly fast decline in velocity dispersion. 
While the prospects for spatial resolution in 
the dispersion are not good, it can be used to 
constrain possible models for the galaxy: an 
isothermal mass distribution looks problem¬ 
atic. The models of Keeton and Kochanek 
(1997) explore some of these possibilities. A 
hnal possibility is that the shear from cluster 
is pushing the images closer together, over¬ 
coming some of the splitting from the lens 
itself, but it is unclear whether such a model 
could be workable. 

In recent years the sophistication of grav¬ 
itational lens models has clearly outstripped 


the quality of the data available. The Keck 
telescopes and the LRIS spectrograph are 
truly marvelous facilities, and it is the aim 
of this and subsequent papers to try to bring 
some grist to the hungry mill, and perhaps 
help teach us about the cosmology of our uni¬ 
verse. 

Thanks are due to Paul Schechter for en¬ 
couragement and helpful suggestions through¬ 
out this project. 



Table 1 

Observing Log. 


Obs# 

Objects 

UT Date 

UT 

sec 2 : 

PA 

Exposure 

Slit 

Grating 

1. 

AGK2+14783 

3/30 

5:46 

1.01 

90 


0.7 

600/5000 

2. 

PG1115 GL+B 

3/30 

8:55 

1.02 

35 

1500 

1.0 

600/5000 

3. 

PG1115 GL+B 

3/30 

9:22 

1.02 

35 

1500 

1.0 

600/5000 

4. 

PG1115 Al+Al 

3/30 

9:53 

1.04 

35 

500 

1.0 

600/5000 

5. 

PG1115 G1+G2 

3/30 

10:10 

1.05 

42 

1500 

1.0 

600/5000 

6. 

PG1115 G1+G2 

3/30 

10:36 


42 

1500 

1.0 

600/5000 

7. 

MS1358 G1+... 

3/30 

11:12 

1.37 

24 

1500 

1.0 

600/5000 

8. 

MS1358 G1+... 

3/30 

11:39 

1.36 

24 

1500 

1.0 

600/5000 

9. 

B1422 G2+G3 

3/30 

12:21 

1.00 

67 

1500 

1.0 

300/5000 

10. 

B1422 G2+G3 

3/30 

12:47 


67 

1500 

1.0 

300/5000 

11. 

B1422 GL 

3/30 

13:54 

1.09 

33 

1500 

1.0 

300/5000 

12. 

B1422 GL 

3/30 

14:21 

1.15 

33 

1500 

1.0 

300/5000 

13. 

HD132737 

3/30 

15:11 

1.20 

90 


1.0 

300/5000 

14. 

HD132737 

3/30 

15:23 


90 


0.7 

600/5000 

15. 

AGK2+14783 

3/31 

5:26 

1.01 

90 


0.7 

600/5000 

16. 

AGK2+14783 

3/31 

5:45 

1.01 

90 


1.0 

300/5000 

17. 

PG1115 GL+B 

3/31 

8:11 

1.04 

35 

1500 

0.7 

300/5000 

18. 

PG1115 GL+B 

3/31 

8:35 

1.02 

35 

1500 

0.7 

300/5000 

19. 

PG1115 A1+A2 

3/31 

9:05 

1.02 

35 

500 

1.0 

300/5000 


Note.— PA is east from north, exposures are in seconds, slit widths are in arcseconds. 
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Table 2 

Redshifts and Dispersions. 


Galaxy 

y 

z 

± 

a 

± 

r 

B1422-GL 


0.3366 

0.0004 



4.2 

B1422-G2 

-2.9 

0.3376 

0.0001 



13.0 

B1422-G3 

+2.9 

0.3427 

0.0001 



11.4 

B1422-GX 

-1.9 

0.5360 

0.0005 



em. 

PG1115-GL 


0.3098 

0.0002 

281 

25 

10.4 

PG1115-G1 

-6.7 

0.3098 

0.0001 

256 

20 

14.8 

PG1115-G2 

+6.7 

0.3123 

0.0001 

130 

60 

10.9 

PGlllS-Gx 

+2.9 

0.3121 

0.0002 



em. 

PG1115-GL 


0.3095 

0.0002 



8.1 

MS1358-G1 

+9.6 

0.3272 

0.0001 

299 

22 

11.9 

MS1358-G2 

-9.5 

0.3235 

0.0001 

235 

23 

14.1 

MS1358-G3 

-1.7 

0.3248 

0.0001 

155 

56 

9.4 

MS1358-G4 

+18.3 

0.3229 

0.0002 

175 

86 

5.2 


Note.— Columns: Galaxy name, slit position ("), red- 
shift and error, velocity dispersion (kms“^) and error, 
cross-correlation r value. 
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